NASA Technical Paper 1754 


Surface Pyrometry in Presence 
of Radiation From Other Sources 
With Application to /Turbine- 
Blade Temperature Measurement 

' ' } 

■■ ■ \ ' ,4 

Donald R, Buchele 

NOVEMBER 1980 

■ I 


NASA 
TP 
1754 
c. 1 


I 


h' 


loan con 
afwltech 

mi/iio A 


as 

b-" ? 

!^i 

SI 

LJ! 




TECH LIBRARY KAFB. NM 



NASA Technical Paper 1754 


tech library kafb, nm 



□13MaS3 


Surface Pyrometry in Presence 
of Radiation From Other Sources 
With Application to Turbine- 
Blade Temperature Measurement 


Donald R. Buchele 
Lewis Research Center 
Cleveland, Ohio 


IWNSA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1980 



Summary 

Surface pyrometry can be used to measure surface- 
temperature distributions on rotating turbine blades 
in jet engines. However, advanced turbofan engines, 
operating at high gas temperatures and pressures, can 
produce interfering gas radiation greater than the 
radiation from turbine blades. To measure and cor- 
rect for interfering radiation, several methods were 
compared that use multiple-wavelength pyrometry. 
The methods included blade temperature reduction 
by increased blade cooling and use of one reference 
thermocouple on the blade. It was concluded from an 
error analysis of each method that measurement at 
two wavelengths was best for blade temperature 
determination. Accuracy of the measured 
temperature could be substantially improved by an 
auxiliary measurement, at the same two wavelengths, 
in which the pyrometer was sighted on a portion of 
the blade to which a thermocouple was attached. The 
conclusions are of general applicability to surface 
temperature measurement and not Just to turbine- 
blade temperature measurement. 


Introduction 

Surface-temperature determination by measure- 
ment of the radiant energy emitted from the surface 
may be subject to errors caused by the presence of in- 
terfering radiation. A typical application in which 
this problem arises is the measurement of 
temperature distribution over the surface of a cooled 
turbine blade in a jet engine (ref. 1). Although the 
methods treated in this report are applicable to the 
general problem of surface-temperature measure- 
ment, the treatment will be directed specifically to the 
measurement of turbine-blade temperature because 
this measurement is of immediate concern (refs. 2 
and 3) and because such specificity permits a quan- 
titative estimate of attainable accuracy. 

Tests have shown that thermal radiation emitted 
by turbine blades is augmented by radiation from 
gas, flame, or particulates in the line of sight of a 
pyrometer or by reflection of this radiation by the 
blades into the line of sight. This interfering radia- 
tion causes an error in blade temperature measure- 
ment. The magnitude of the error depends on blade 
temperature, blade emittance, gas temperature, and 
gas emittance. Correction for interfering radiation in 


other applications is considered in references 4 to 7 . 
In application to a jet engine (as reported in ref. 1) 
radiance measurements at two wavelengths were used 
to determine blade temperature and effective gas 
emittance, with independently determined blade 
emittance and gas temperature. 

In this report methods are considered that depend 
less on independently determined parameters. An er- 
ror analysis of the methods enables their comparison. 
The results of the comparison provide a selection of 
methods that should lead to minimal error in blade 
temperature determination. 


Analysis 

Definition of Radiation Parameters 

A pyrometer, sighted on a turbine blade in an 
engine, can receive blade radiation plus interfering 
radiation that is reflected by the blade from other 
sources. These sources include the surrounding 
blades, gas, flame, or particulates. The gas and the 
particulates in it are usually at a much higher 
temperature than is the blade itself. The radiation go- 
ing from the blade to the pyrometer can be partially 
absorbed by the gas and by the particulates and sup- 
plemented by emission from the gas and the par- 
ticulates. Considering all these sources of radiation, 
reflection, and absorption, it is reasonable to 
separate the observed radiation into two com- 
ponents: one from a source at a relatively low 
temperature that represents the blade radiation, and 
one from a source at a relatively high temperature 
that represents the effect of the gas and particulate 
radiation. The radiation originating from the blades 
can be represented by two parameters: a temperature 
Tb, and an effective emittance ex,6- The radiation 
originating from the gas (including flame and par- 
ticulates) can be represented by two parameters, an 
effective temperature Tg and an effective emittance 
e\,g. (All symbols are defined in appendix A.) The 
radiation received by the pyrometer has a spectral ra- 
diance L\ given by the sum of the two component 
radiances 


-t- Tx,g (i) 



For error analysis the Planck radiation equation can 
be replaced by the Wien approximation. When the 
Wien equation is used, equation (1) becomes 

= CiX-5 [ + ex.ge-"2/X7>] (2) 

The quantity is an effective emittance that in- 
cludes the effect of radiation from adjacent blades 
that is reflected from the blade under observation. 
The quantities Tg and e\^g are an effective 
temperature and an effective emittance, respectively, 
that include the effects of gas or particulate radiation 
that is in the line of sight or that is reflected from the 
blade into the line of sight. 

Radiation from the water vapor component of the 
gas and radiation from the turbine blade are com- 
pared in table I for a turbojet and for an advanced 
turbofan engine. Engine operating conditions are 
listed at the combustor, upstream and downstream of 
the first row of turbine blades, and at the turbine 
blade. The emittance of water vapor at 1 nm 
wavelength is calculated by using the gas absorption 
coefficient from reference 8. The product of e\ and 
the blackbody radiance L\ at the listed temperature 
gives the radiance e\L\. The last column in table I is a 


ratio of gas radiance to blade radiance. With a 
10-centimeter path length the gas radiance is com- 
parable to the blade radiance for an advanced tur- 
bofan. Blackbody radiance is very large at the com- 
bustor temperature. Combustor radiation is shown in 
reference 1 to propagate downstream past the turbine 
inlet guide vanes and then to be reflected from the 
turbine blade at a level comparable to that of the 
blade radiation itself. 

The computation of infrared radiant emission and 
absorption from 1 to 10 fim by combustion gases is 
extensively treated in reference 9. Because there is 
relatively little radiant emission from 0.6 to 1 /xm, 
there is also a lack of published data on gas absorp- 
tion constants in this band. Although the weakness 
of radiation in this region is an advantage for 
turbine-blade pyrometry, the lack of data prevents 
estimation of the gas spectral emittance. Actual 
measurement of the spectral emittance is required for 
each experimental situation. This stringent require- 
ment can be simplified somewhat by eliminating the 
need to determine gas temperature and emittance as 
separate quantities. The two variables ex,g and Tg in 
equation (2) can be replaced by Lx,g from equation 
(1) to give 


TABLE I. - OPERATING CONDITIONS IN A TURBOJET AND AN ADVANCED TURBOFAN ENGINE 


[Gas emittance estimated for HgO at 1 fxmj HgO partial pressure, 0.13 of total pressure; viewing 
path length, 1 0 cm . ] 



Static 

temperature, 

K 

Total 

pressure, 

N/m^ 

Emittance, 

Blackbody 

radiance, 

9 

W/cm cm sr 

Radiance, 

^X^X 

Ratio of gas 
radiance to 
blade 
radiance, 

(^x^x)b 

Turbojet; 







Combustor gas 

2150 

4.7X10^ 


153 000 



Gas upstream of 

1150 

4.7 

0.0018 

450 

0.83 

0.003 

turbine blade 







Gas downstream 

920 

1.7 

.00077 

22 

.012 

0 

of turbine blade 







Turbine blade 

1100 


1.0 

260 

260 


Advanced turbofan: 







Combustor gas 

2590 

37 


460 000 



Gas upstream of 

2200 

20 

.0053 

177 000 

940 

3.6 

turbine blade 







Gas downstream 

1810 

17 

.0045 

43 000 

193 

.7 

of turbine blade 







Turbine blade 

1100 


1.0 

260 

260 
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For combustion gases where the emittance is much 
less than unity, gas radiances at two wavelengths, if 
once determined during a test, usually maintain the 
same radiance ratio L 2 ,g/L\^g to a higher degree than 
the absolute radiance at any one wavelength. For 
solid surfaces like those of turbine blades, the ratio 
of monochromatic emittances at two wavelengths 
and the same temperature are known or calculable 
with a higher percentage of accuracy than the ab- 
solute value of emittance itself. (When the emittance 
approaches unity as the upper limit because of the 
cavity effect of adjacent turbine blades, a ratio of 
reflectances at two wavelengths, that is, 
(1 -e 2 , 6 )/(l -ei,*), may be used instead of a ratio of 
emittances.) Consequently, where several equations 
of the form of equation (3) are to be solved 
simultaneously, the accuracy of the calculations is 
improved if the ratios are used in place of absolute 
values, as far as possible. This approach is tanta- 
mount to recognition that the quantities L\^g, Z. 2 ,g, 
and Z, 3 _g are not truly independent variables but are 
related to each other. Similarly, the quantities 
e 2 ,b, and e-i,b are not truly independent variables but 
are related to each other. The ratios r,,* and can 
be defined by the equations 

r2,b-^2.b/^\,b 

•^3,g=^3,g/Ti_g 

Equation (3) written for three wavelengths thus 
becomes 


Z„=CiXi5e,^-^/Xl7A + Z.,,g 
L3 = 

Radiance is measured with an instrument having a 
bandpass spectral response. The error analysis 
assumes a sufficiently narrow bandpass for equations 
(4) to be defined for a monochromatic wavelength 
near the center of each pass band. 

Table II lists values of the blade emittance ratios 
actually measured for some used blades taken from 
test engines. It must be emphasized that the ratios of 
effective emittances that appear in equations (4) are 
affected by blade surface deposits and the cavity 
effect of adjacent blades when blades are installed in 
an engine. These ratios are necessarily less well 
known than the ratios listed in table II, since these 
latter ratios were measured when the blade was 
isolated in a laboratory environment. 

Radiometric Methods 

In this discussion, we distinguish between the 
direct radiometric measurement obtained with the 
pyrometer and the independent determination of 
various parameters. The latter can be obtained by 
analytic computation, by a separate laboratory ex- 
periment, or by a measurement in situ with auxiliary 
instrumentation. 

With a measurement of L\ at one wavelength, 
equation (3) can be used to determine 7), when 
independent measurements are made of e\,6 and 
Lx,g. The total error of Tb is caused by the error in 
radiometric measurements of L\ and the error of 
each independent measurement. By virtue of equa- 
tion (2) it is suggested that, when gas and blade 



T.ABLE II. - EXAMPLES OF EMITTANCES AND EMITTANCE RATIOS 
[Used blades measured at NASA, at 300 K blade temperature.] 


Blade 

Wavelength, 

K, pm 

Emittance ratio 


Xg = 0.75 

X 3 = 0.55 

^ 2 ,b/^l,b 

^3,b/^l,b 


Emittance, 


B-1900 

0,855 

0.875 

0.905 

1.02 

1.06 

B-1900, later test 

.91 

.92 

.93 

1.01 

1.02 

Blade coated with 

.53 

. 555 

.58 

1.05 

1.09 

ZrSiO^ certtmlc 
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temperatures are different, measurement at two 
wavelengths would be effective in discriminating be- 
tween blade and gas radiances in equation (1). Fig- 
ure 1 shows blade and gas radiances in a spectral 
range from 0.55 to 1.0 ixm, with typical temperatures 
7* of 1 100 K and Tg of 2000 K, The ratio L\^g/L\^t 
of radiances substantially changes with wavelength. 
The blade radiance is a larger fraction of the total ra- 
diance at the longest wavelength; this wavelength is 
therefore most effective for determining blade emit- 
tance and temperature. At the shorter wavelength the 
gas radiance is a larger fraction of the total radiance; 
this wavelength is thus most effective in determining 
gas radiance. These conclusions about the relative 
effectiveness of the shorter and longer wavelengths 
assume that the gas spectral emittance is a constant 
value, such as that produced by particulates at the 
gas temperature. If the gas spectral emittance were to 
become 1/7 as great at 0.75 (im as it is at 1 /tm (which 
is possible with gas molecular radiation), the equa- 
tions at two wavelengths would become redundant, 
and their simultaneous solution would be indeter- 
minate. Thus two wavelengths of measurement are 



not usable when the ratio L\^g/L\^b of radiances is 
the same at both wavelengths. 

The methods considered herein are those that may 
reduce the number of independently determined 
parameters or the degree of dependence of the 
temperature determination on those parameters. The 
methods are combinations of four different measure- 
ment conditions: 

(1) Measurements at several wavelengths, A — Ra- 
diance measurements are made at two or three 
wavelengths. This method uses two independently 
determined parameters ex,f, and L\^g at each 
wavelength, but only one radiance measurement is 
made at each wavelength. 

(2) Measurements at two blade temperatures, 
B — increased blade cooling is used to produce a 
lower blade temperature with the same gas 
temperature. This method introduces one additional 
parameter, the second blade temperature Ty and an 
added radiance measurement L\ at each wavelength. 

(3) Thermocouple on the blade, C — A thermocou- 
ple is used to measure temperature at one point on 
the blade. The additional parameters L\^g and 
found when viewing the thermocouple junction can 
then be applied to determine blade temperature at 
other locations on the blade. The error in measuring 
temperature with the thermocouple is assumed to be 
negligible in comparison with the errors of optical 
pyrometry as here used. 

(4) View a nonradiating surface, D — A low- 
temperature surface or a porthole is used to provide a 
surface with negligible radiation. Only gas radiation 
is measured. Since the line of sight does not intersect 
a blade, the gas radiation parameters may become 
different when the line of sight is directed at the blade 
during a test. 

Some combinations of measurement conditions 
are listed in table III as methods 1 to 8. The column 
“Radiance measurements” lists the number of 
wavelengths and temperatures at which radiance 
measurements are made. The number of equations 
limits the number of parameters that can be solved 
for. The remaining parameters must be determined 
by independent, separately made measurements or 
estimates. When the number of equations exceeds the 
number of parameters to be solved for, a least- 
squares method of solution is to be used. 

The methods are conveniently arranged in two 
groups. Primary methods are those that alone would 
suffice to determine 7), if there were adequate 
knowledge of the listed “Parameters to be in- 
dependently determined.” Auxiliary methods are 
those that are used in conjunction with one of the 
primary methods to provide an actual experimental 
determination of one or more of the parameters that 
must be independently determined, with resultant im- 
provement in accuracy. The result of such combina- 
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TABLE m. - RADIOMETRIC METHODS 


(a) Primary methods for determining blade temperature Tj^ 


Method 

Measure- 

ment 

condi- 

tions 

Radiance measurements 

Number 

of 

Parameters solved for. 

Parameters 

independently 

determined, 

Qk 

Wavelength, 

Temperature, 

K 

tions 

Blade 

Gas 

Blade 

Gas 

1 

A 




1 

^b 


^l.b 

^l.g 

2 

A 



^b 

2 

Tb 

^l.g 

^l,b’ ^2,b 

^2,g 

3 

A 



Tb 

3 

Tb 

^l.g’®2,g 

^l,b’’’2,b’’^3,b 

^3,g 







♦ 




4a 

A, 

B 


■rb-Tb 

4 


^l.g 

^l,b’ ’'2,b 

«2,g 







♦ 




4b 



Xj^.Xg 

Tb’Tb 

4 

Tb-Tb 

^l,g*^2,g 

^l,b’’^2,b 

— 

5a 



Xj^, Xg, Xg 

Tb 

5 

Tb>< 

^l.g’ ^2,g 

®l,b’’^2,b>’^3,b 

«3,g 




Xi.Xg 

< 






5b 




Tb 

5 


^,g’ ^2,g’’^3,g 

^l,b’’’2,b’’^3,b 

— 




Xj^.Xg 

< 







(b) Auxiliary methods 


6 






^X,g’^X,g 
^X,g’ ^X,g 



A,B, C 


* 

2 


%b 

T..T.* 




^b’ ^b 


^X,b’ ’’x.b 


8 

A,D 





^X,g’^X,g 










(c) Combinations of methods 


2&6 

A, C 


Tb 

4 

Tb 

^l.g’ ^2,g 

^l,b’ ’'2,b 
♦ 

— 

2&7 

A,B, C 

^1’ ^2 

Tb-^b 

6 

''’b’^l,b* ’^Z.b 

^l.g’«2,g 


— 

2&8 

A,D 

^l’^2 

Tb 

4 

Tb 

^i,r ^2,g 

®l,b’^2,b 

— 


lions with primary method 2 is listed in table III(c). 
Method 2 was selected as the best primary method. In 
some combinations of primary and auxiliary 
methods accuracy is further improved if the radiance 
measurements required for the auxiliary method are 
the same radiance measurements that were used in 
the primary method. 

The methods of solution for Ti, are 

(1) Measurement at one wavelength, using one of 
equations (4) — The interfering radiation Z.i,g must be 
determined by an auxiliary method. 

(2) Measurement at two wavelengths, using two of 
equations (4) 

(3) Measurement at three wavelengths, using three 
of equations (4) 

(4&5) Methods 2 and 3 repeated with additional 
measurements at a second, lower blade 
temperature — The methods differ in the number of 


parameters analytically determined. Methods 4b and 
5b solve for all gas radiation parameters. 

The auxiliary methods listed in table Ill(b) are one- 
wavelength measurements repeated at other 
wavelengths as needed. The methods are 

(6) Measurement with a thermocouple on the 
blade, giving T, at that point — This is an inverse of 
method 1, where Tt, is solved for. When this method 
is applied at two or three wavelengths, the quantities 
L\^g are expressed as ratios /? 2 ,g and Ri,g. Thus all in- 
dependent gas parameters needed by methods 1 to 5 
are determined. Some gain in accuracy may be ex- 
pected by operating at minimum blade temperature 
in order to minimize blade radiation. 

(7) Measurement with a thermocouple on the blade 
and with blade cooling to give two blade 
temperatures T, and T ' — The use of two blade 
temperatures allows determination of blade emit- 





tance; method 6 does not. Thus all independent 
parameters needed by methods 1 to 5 are determined. 
Where the conditions of the experiment permit, some 
gain in accuracy can be expected by applying this 
method in two steps. By operating at a test condition 
with small gas radiation, the blade emittance is more 
accurately determined. Then, by operating at the 
normal test condition where gas radiation is larger, 
L\g is more accurately determined. 

(8) Measurement along a line of sight to a surface 
of negligible thermal radiation and negligible reflec- 
tion of incident radiation — This method determines 
gas radiation along the line of sight, but it neglects 
gas radiation that may be reflected from the blade 
when the blade is viewed by methods 1 to 6. Thus the 
ratios /? 2 ,g and may be inaccurate when reflected 
gas radiation is significant. 

The combinations of methods listed in table III(c), 
when compared with the primary method 2, show a 
transfer of some parameters from the column 
“Parameters independently determined” to the col- 
umn “Parameters solved for.” The resulting im- 
provement in the accuracy of determining Ti, is later 
derived by the error analysis and calculated for a 
specific example. 

Error Estimate for Radiation Parameters 

In the following discussion, an error of known 
algebraic sign in a variable x is denoted by dx, and a 
random error in x is denoted by 5x. 

The parameters solved for in table III have errors 
caused by 

(1) The error of radiance measurement at each 
wavelength (type 1 error) 

(2) The error of the parameters independently 
determined in table III (type 2 error) 

In order to compute the errors by using matrices, 
let Pj represent parameters solved for, let Qk repre- 
sent parameters to be independently determined, and 
let Lj represent measured radiances at wavelengths X,-. 

Type 1 errors are derivable from 

where (5) 

The quantity IP/y is called the error sensitivity coeffi- 
cient. The nondimensional partial derivatives are ob- 
tained by inverting the matrix of nondimensional 
partial derivatives (3I,;/P/)(P//L/). These derivatives 


are found with a computer program by numerically 
differentiating each equation for with respect to 
the parameters Pj listed in table III for each method. 

Methods 4a and 5a in table III have more equa- 
tions than unknowns. The excess data are accom- 
modated by an additional matrix manipulation that 
uses the least-squares method (ref. 10). 

Type 2 errors are derivable from 






The quantity Xj^k can be expressed differently by 
noting that 



( 7 ) 


where, under the summation sign, the first partial 
derivative is obtained from the solution of the 
simultaneous equations (4) but the second partial 
derivative is obtained from that one of equations (4) 
which represents £,■ so that 

/= 1 


Many of the terms in equation (8) are zero in practice 
since, for example, the partial derivative of the ra- 
diance at the i'h wavelength with respect to the emit- 
tance at the kth wavelength (i.e., dLj/dQi^, where 
Qk = ^k,b) is zero. 

The expression for the total random error in Py due 
to random errors in both L,- and Q/c depends on 
whether any of the Qk(fc=l to /, where l<m) are 
determined by an auxiliary method that uses the same 
radiance measurements as the primary method. For 
such Qk we note that 



where 


(Ar=l,2,...,0 


( 9 ) 


y _ ^Qk Lj 

^Li Qk 


( 10 ) 


so that the total error is 
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Then the total variance in Pj is given by 



where W^ij, and Yi^k are given by equations (5), 
(8), and (10), respectively. 

Numerical Solutions and Selection of 
Methods 

The nondimensional differential coefficients in 
equations (5) and (6) will be determined for a 
representative example where sources of error are 
sufficiently large to provide a clear indication of the 
relative merits of the various methods in securing ac- 
curate determination of 7*. 

Assumed temperatures are 1 100 K for the normally 
cooled blade, 1000 K for the blade when increased 
cooling is applied, and 2000 K for the gas. Assumed 
center wavelengths of measurement are 1, 0.75, and 
0.55 fim; these are obtainable and measurable with 
bandpass filters and a silicon radiation detector. 
Subscript 1 applies to the longest wavelength, 
subscript 3 to the shortest. 


Solutions for Tb by Methods 1 to 5 

The error sensitivity coefficients for blade 
temperature are shown in figure 2 and table IV. In 
figure 2 the coefficient (eq. (5)) is caused by radiance 
measurement error; in table IV the coefficient (eq. 
(6)) is caused by an error of independently deter- 
mined parameters. The abscissa L\^g/L\^b bi figure 2 
is a ratio of gas radiance to blade radiance at 1 /rm. 
The ratio includes the effects of blade and gas emit- 
tances. In table IV the error sensitivity coefficient for 
blade emittance is nearly independent of L\^g/L\^b> 
but the error sensitivity coefficient for gas radiance is 
nearly proportional to L\^g/L\^b- The total blade- 
temperature error indicated by each method is given 
by the square root of the sum of the squares of errors 
computed with the error sensitivity coefficients in 
figure 2 and table IV. 

Solution for t\b ^nd L\_g by Methods 6 to 8 

The error sensitivity coefficients for ex,* and 
are shown in figure 3 and table V, which parallel the 
previous figure 2 and table IV in definition of terms. 
As these methods are one-wavelength methods, the 
results are shown for each of the three wavelengths. 
In figure 3 the coefficient (eq. (5)) results from a ra- 
diance measurement error; in table V the coefficient 
(eq. (6)) results from an error in blade temperatures 
Tt and T* as derived from thermocouple readings and 

Method 



lljm, 


.0001 .001 .01 
Gas emittance at 1 gm, g 

Figure 2. - Error sensitivity coefficients for 
methods 1 to 5. Temperature of biade with 
increased cooiing 1000 K; temperature 
of cooied blade, Tj,, 1100 K; gas temperature. 
2000 K; blade emittance, ex_[,, 1. 
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TABLE V. - AUXILIARY METHODS, ERROR SENSITIVITY COEFFICIENTS CAUSED BY 


INDEPENDENTLY DETERMINED PARAMETERS WHEN L, „/L, , = 1 

[When dependent variable is gas radiance, the coefficient is to be divided by the actual value of 
Lj^ g /^1 b‘ dependent variable is emittance, the coefficient is independent of 


Method 

Wavelength, 

X, 

gm 

Error sensitivity coefficient 


* + 

9^X.b 

^^X,b,., '^t 


3T. L, „ 

t X, g 

®^t ^X,g 

®^X,b ^X,g 

®"^t ^X,b 

^■^t ^X,b 

6 

1.0 

-13.0 

(a) 

-0.99 

(a) 

(a) 


.75 

-2.43 

(a) 

-.14 

(a) 

(a) 


.55 

-.19 

(a) 

-.01 

(a) 

(a) 

7 

1.0 

4.8 

-5.3 

(a) 

-17.9 

5.3 


.75 

.5 

-.6 

(a) 

-21.1 

4.1 


.55 

.02 

-.02 

(a) 

-26.2 

2.7 


Quantity does not exist. 


dents for ei^t, are of similar magnitude for all 
methods listed. Thus methods 2 to 5 are equally 
affected by the error of the independently determined 
blade emittance. 

In the absence of other strong reasons for selection 
among the methods, selection may be made on the 
basis of minimizing the number of measurements. 
Since method 2 requires the fewest number, it was 
selected. 

Auxiliary methods.— li now remains to choose 
among method 2 alone and combinations of method 
2 with one of the auxiliary methods 6, 7, and 8. 
Table VI lists the values of the error ratios that enter 
into the use of each of these method combinations, 
for four values of the ratio L\^g/L\^i,. Table VI lists, 
in the second column, the square root of the sum of 
the squares of all the error coefficients that have ra- 
diance as an independent variable; this number 
represents the total error, in percent, caused by a 
1 percent error in each radiance measurement. The 
remaining columns list the error coefficients for in- 
dependently determined emittance parameters. 

If the combination of methods 2 and 6 is used, the 
radiance measurements used in method 6 should be 
the same radiance measurements that were used in 
method 2. If this is done, errors due to an incorrect 
assumption of ex,* cancel out; hence, an accurate 
determination of ex,* is not necessary for the portion 
of the blade near the thermocouple. This statement is 
no longer fully accurate for regions of the blade 
where blade emittance is distinctly different from the 
emittance at the thermocouple. In such an event the 


cancellation of errors is not complete and a 
knowledge of the change in 6x,* is desirable to im- 
prove accuracy. This situation is not the situation 
treated in table VI. 

Those combinations of methods where the variable 
under consideration is one of the variables solved for 
are indicated by footnote (a) in table VI. 

When method 6 or 7 is used, any error in blade 
surface-temperature measurements with the ther- 
mocouple will introduce an equal error in blade 
temperature measurement. When method 8 is used, 
no error cancellation similar to that in method 6 
occurs. However, because the gas radiance is 
measured, the error ratio dTb/dR-i,g is omitted, and 
the effects of radiance errors in determining are 
included in the error ratio dTi,/dL. 

The coefficient values in table VI were calculated 
with the same ratio Z.i,g/Z.;,* for both method 2 and 
the auxiliary methods. When these methods are ap- 
plied sequentially, with different engine conditions, 
the ratio Z.i,g/Z,i,* will be different for each method. 
This will change the values of the coefficients listed in 
table VI for the auxiliary methods. It is evident from 
figure 3 that smaller values of the coefficients will be 
obtained if ex,* is determined by method 7 when 
L\,g/L\^b is small and if L\g is determined by 
methods 6 or 7 when Z.i,g/Li,* is large. Whenever 
possible this procedure should therefore be used. 

Method 7 can be used for an initial determination 
of blade emittance, which is then used with method 6 
for the majority of data taking. This procedure will 
reduce the time spent in cooling the blade. 
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TiiBLE VI. - TOTAL ERROR RATIOS FOR BLADE TEMPERATURE, WHEN 


DATA FOR METHOD 2 ARE SUPPLIED BY AN AUXILIARY METHOD 


Radiance 

Total error ratio 

ratio, 

^l,b 



^2 e 

dTb 

„ ^l.b 

dTb 

,.*’2,b 

dTb 

^t 


‘^«2,g 

Tb 

‘^^l.b 

" Tb 

^^2,b 

" Tb 

dTt 

Tb 

Method 2 

0.01 

0.096 

0.0009 

-0. 

094 

0.0132 


.1 

.106 

.0094 





■ 

■ 

1 

.217 

.0948 





■ 


10 

1.42 

.9480 





■ 

■ 

Methods 2 and 6*^ 

0.01 

0.135 

(a) 

0 

0 

■ 

■ 

.1 

.150 







■ 


1 

.307 







H 


10 

2.01 







■ 

■ 

Methods 2 and 7^ 

0.01 

0.135 

(a) 

(a) 

(a) 

1 

. 1 

.176 









1 

10 

.354 

2.29 



1 





f 

' 


Methods 2 and 8 

0.01 

0.096 

(a) 

-0. 

094 

0.0132 

n 

.1 

.106 







■ 


1 

10 

.255 

1.95 







1 

1 


^Quantity does not exist. 

Entries apply only where blade emittance is approximately equal to blade 
emittance near the thermocouple. 


The accuracy of method 7 depends on the 
temperature difference T* - produced by blade 
cooling. The effect of this temperature difference on 
the values of blade emittance error coefficients is 
shown in table VII for two values of the ratio 
L\^g/L\^b and for Xj, X 2 . Large values of the coeffi- 
cients are found for blade emittance at X 2 when 
L\,g/L\,b=l. Small values are found when gas ra- 
diance is small, at L\^g/L\^b = 0. 1 . A temperature dif- 
ference as small as 25 K can give acceptable accuracy. 


The principal interest of this report has been 
maintenance of pyrometer accuracy in the presence 
of interfering gas radiation. This accuracy is always 
reduced by the presence of interfering radiation. It is 
thus important to demonstrate acceptable accuracy 
of the pyrometer under the ideal condition where 
there is no gas radiation. Accuracy of calibration and 
measurement is affected by the wavelength assigned 
to the bandpass filters. These topics are treated in 
appendix B. 
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TABLE VII. - ERROR SENSITIVITY COEFFICIENTS FOR AUXILIARY 

METHOD 7 


[ Coefficients are rms of coefficients at T^. and T^^ . ] 


Error sensitivity 

Radiance 

Wave- 

Temperature difference between 

coefficients 

ratio. 

length. 

turbine blade and blade with 



X 

increased cooling 
* 

» 


^l.b 



Tb-^b- 

K 





25 

50 

100 

200 




0.1 

1 

5.3 

2.7 

1.6 

1.1 




.1 

2 

6.6 

3.7 

2.3 

1.9 


^X,b 


1 

1 

59 

33.1 

20.7 

15.8 


1 

2 

307 

181 

123 

104 




0.1 

1 

62.4 

32.3 

18.7 

13.9 




.1 

2 

63.3 

33.9 

21.5 

17.8 

1 

*^^t ^X,b 

1 

1 

1 

62.4 

32.3 

18.7 

13.9 


2 

2 

63.3 

33.9 

21.5 

17.8 


Error Estimate for Blade Temperature 

The error sensitivity coefficients that were 
presented in table VI for several methods are con- 
verted to blade temperature error by using assumed 
errors for radiance measurement and for blade and 
gas emittance parameters. The blade temperature 
error given by these methods is compared with the 
error of a one-wavelength pyrometer when gas radia- 
tion is ignored. (This is method 1 with an error 

bLi\^g = 

The blade temperature error is listed in table Vlll 
for each assumed error, for four values of the ratio 
L\_g/L\^f}- Table VIII lists, in the second column, the 
error caused by radiance measurement and, in the 
next four columns, the error caused by incorrectly 
estimating blade emittance and gas radiance. The last 
column is an rms total error of the preceding col- 
umns. The assumed errors listed in the column 
headings are obtained as follows: 

(1) Radiance measurement error is bL}^/L\ = 0.Q2 

(2) From table II, blade emittance e\^i, = 0.9l, with 
error estimate 6ei,b = 0-05, gives 5ei,6/ei^6 = 0.055. 

(3) From table II, the ratio ri,b = with error 
estimate 5r2,6 = 0.01, gives 5r2,*/r2,* = 0.01. 

(4) For method 1 only, where L\^g is ignored, 

bL\^g = Li^g. 

(5) The ratio /?2,g depends on the emitter. For par- 
ticulate matter the emittances at two wavelengths are 
nearly alike. For molecular emitters the emittances at 
two wavelengths may be substantially different. The 


ratio also depends on the gas temperature through 
the Wien equation, as shown by equation (1) and as 
plotted in figure 1 for the 2000 K gas. Because of the 
considerable variation in the ratio with engine 
operating conditions, an error 6f?2,g/i?2,g = 0.3 is 
assumed. 

Comparison of the rms total error for each method 
shows that method 2 is not better than method 1, for 
the assumed error 6i?2,g/^2,g- To show some im- 
provement, it is necessary to determine /?2,g with less 
than 30 percent error. The auxiliary methods are 
effective in doing this. Methods 2 and 6, and 2 and 7, 
also reduce the blade temperature error caused by 
blade emittance error. These methods also improve 
the accuracy when gas radiation is small. 


Summary of Results 

Several methods of using a multiple-wavelength 
pyrometer to measure interfering radiation have been 
shown by an error analysis to permit improvement in 
the accuracy of blade surface-temperature measure- 
ment. It was concluded that measurement at two 
wavelengths was best for blade temperature calcula- 
tion. Accuracy of the calculated temperature was 
substantially improved by an auxiliary measurement 
at two wavelengths with the optical line of sight on an 
element of the blade surface where a thermocouple 
was attached to the blade. Blade emittance can also 
be determined by repeating the auxiliary 


II 





TABLE Vm. - COMPARISON OF BLADE TEMPERATURE ERROR BY SEVERAL METHODS 
USING ESTIMATED ERRORS OF RADIATION PARAMETERS 



Quantity does not exist. 

measurements after cooling the blade to obtain a Lewis Research Center, 

temperature drop greater than 25 K. The methods National Aeronautics and Space Administration, 

discussed herein can be generalized to other applica- Cleveland, Ohio April 15, 1980. 

tions of surface-temperature measurement in the 505-04. 

presence of interfering radiation from other sources. 
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Appendix A 
Symbols 


Cl 

Planck’s radiation constant, 3.7403 x 10-12 
W cm2 

X 

Y 

error sensitivity coefficient, eq. (7) 
error sensitivity coefficient, eq. (10) 

C2 

Planck’s radiation constant, 1.4388 cm K 

e 

emittance 

K 

slope of line on graph 

X 

wavelength 

L 

P 

radiance 

primary parameter to be solved for 

T relative transmission factor of filter 

Subscripts: 
b turbine blade 

Q 

R 

independently determined parameter 
ratio of gas radiance at two wavelengths 

r 

ratio of blade emittance at two wavelengths 

s 


S 

relative spectral response of detector 

T 

Ti 

temperature 
calibration temperature 

T 

temperature 

t 

thermocouple 

Ti 

V 

calibration temperature 
measured voltage at temperature T 

1,2,3 

wavelengths (1, longest; 3, shortest) 

Vi 

measured voltage at temperature Ti 

Superscript: 

iV 

error sensitivity coefficient, eq. (5) 


turbine blade with increased cooling 
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Appendix B 

Calibration and Computation Procedures 


Calibration by Use of Effective Wavelength 

The radiometer can be calibrated with a blackbody 
at two or more temperatures. The Wien function 
represents a linear relationship between the coor- 
dinates In L\ and l/T when the wavelength is cons- 
tant. Because of the bandwidth of the filter the 
wavelength to be used with the Wien function is not a 
constant but becomes shorter as the temperature 
rises. Thus a calibration curve in the form of a graph 
of In L\ as a function of 1 /Tis almost, but not fully, 
linear. The deviation from linearity will determine 
the number of calibration points required. Since the 
detector output voltage FIs proportional to L\, In V 
is also a nearly linear function of l/T. Thus linear in- 
terpolation can be used to give \/T when In V is 
known. The wavelength to be used with the Wien 
function at any temperature is derivable from the 
slope of the calibration curve and is given by 

1_ 1 rt/lnix] 

An alternative method, using a calibration at only 
one temperature T\ as described in reference II, is 
applied here to illustrate the concept and use of effec- 
tive wavelength. The effective wavelength is an 
average wavelength that replaces the spectral pass- 
band of the pyrometer filter. 

Typical examples of filter transmission and detec- 
tor response are shown in figure 4 . The filter band- 
widths in this example are about 0.1 fim, about 10 
percent of the average wavelength. 



Figure 4. - Filter transmission factors and detector rela- 
tive response. 


The effective wavelength of a filter is dependent on 
temperature. The effective wavelength appears as 
Xj-j-j-when the radiance ratio L\/L\^i or voltage 
ratio V/ V\ is written as 


_F 



(B 2 ) 


where subscript 1 represents values at the calibration 
temperature and Xj-j-j-is the effective wavelength 
for the temperature interval Ti to T. As the effective 
wavelength depends on T\ and T, its defining equa- 
tion will be derived. For a small interval where T=Ti 
the effective wavelength at T, designated Xj-, is a 
weighted average. 


1 _ JJ(£xrS/X)dx 
JJ(TxrS)dX 


(B 3 ) 


where r is the relative filter transmission factor and S 
is the detector relative spectral response. Only Lx 
depends on T. The integration in equation (B 3 ) is re- 
quired for calibration, but it will be avoided for 
temperature computation. A graph of I/X7- as a 
function of l/T is, to a good approximation, a 
straight line that has a slope K between the two 
temperatures T and T\. 


I/X7 — 1 /Xt-| 
l/T-l/Ti 


(B 4 ) 


Over the interval 7 " to T], the average 1/X gives the 
mean effective wavelength (ref. 11) 


1 _ 1 1" 1 1 

^T\-T 2\X7-, \r 


(B 5 ) 


Substituting equation (B 4 ) in (B 5 ) gives the mean 
effective wavelength independent of X7-as 


1 _ 1 K ( I 1 ^ 

Xr,-r Xr/2Vr T, / 


(B6) 
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The slope constant K is calculated by using equation 
(B4) after equation (B3) has been numerically in- 
tegrated by using filter and detector curves in figure 4 
and the Wien function spectral radiance at the two 
temperatures T and Ty. The calibration will provide 
V\ and T\ for equation (B2). 

To calculate an unknown temperature T, given the 
measured voltage V, equation (B2) is solved for Tby 
using \rj _ j given by equation (B6) with an estimate 
of T. By one or two iterations between equations (B2) 
and (B6), the latest value of T is used with equation 
(B6) to get a better value for - r. which gives a 
better value of T" in equation (B2). 

Experimental Check of the Calibration Procedure 

The value of Xj-j based on a calibration with a 
blackbody can be compared with a value calculated 
independently by integrating equation (B3). For ex- 
ample, the effective wavelength of the filters in figure 
4, calculated with equation (B3) at T] = 1 190 K gave 
the following result; 



Filter 

1 

2 

Xj-j, nm 

0.977 

0.756 


An actual radiometer calibration with a blackbody at 
Ti = 1 190 K and one at T 2 = 990 K gave Xj-j _ t with 
equation (B2) and then Xj-j with equation (B6). The 
results are 



Filter 


1 

2 

Vy, volts 

0.467 

0.1629 

V 2 , volts 

0.0381 

0.00652 

\Ty, tim 

^T’l - T’ Mtn 

0.973 

0.757 

0.975 

0.759 


The two methods show good agreement for Xj-j. The 
disagreement for X^j by the two methods is com- 
parable to the difference between Xj-j and Xj-j _ j for 
the temperature interval of 200 K. 

The calibration data were also tested by computing 
the emittance of the blackbody with auxiliary 
method 7 in table III. The result is 



Filter 

1 

2 

‘Kb 

0.988 

0.989 
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